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ABSTRACT 


The consistent problem with operators and planners for power systems has 
been related to minimizing transmission losses. An important role is played 
by reactive power by keeping voltage stability and reliability in the system in 
order to support the transfer of real power. The optimal reactive power 
dispatch is associated with the problem of non-linear optimization along with 
many constraints. In this paper, a study is highlighted for an algorithm that 
optimizes reactive power with the help of particle swarm algorithm and 
compare the result with Newton-Raphson method. Reduction of system 
active power loss is the goal of the function in the projected algorithm. Here, 
the control variables identified are transformer tap positions, generator bus 
voltages, and shunt capacitor banks with switch. This projected algorithm is 
performed on Oman 6 bus electrical grid as oman electricity transmission 
company has an instability voltage issue in chosen 6 bus. 
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1. INTRODUCTION 

Security constraint optimal power flow also abbreviated as SCOPF is demonstrated for obtaining 
reactive power. In the power systems, the operations of the reactive power facilities must be based on the 
certain system-associated technical requirements [1]. These technical requirements are voltage profile 
improvement, reduction in transmission line loss. Also, in order to avoid instability, the voltage stability 
margin is enhanced in order to minimize load perturbation or failure of equipment [2]-[4]. Generally, 
identified issues with power systems that are restructured are procurement of optimal reactive power, 
reduction in losses, and efficient control of voltage. These issues can be expressed as a non-linear 
programming problem (NLP) [5]. The goal is to optimize the formulated function and at the same time 
substantiating the constraints of the objective function. In this objective function, there are two types of 
constraints, equality constraints and Inequality constraints. The equality constraints are generally power flow 
equations whereas the inequality constraints are limits applied on voltages, control variables, line flows, and 
reactive/active power generation [6]. To handle the problems of linear convex, optimization algorithms that 
are based on mathematics are known to be highly efficient. Based on previous studies on reactive power and 
optimized voltage control, the market for reactive power is a non-convex NLP containing multiple feasible 
regions. Hence, an optimum solution can be pointed anywhere within this region of feasibility. Therefore, 
these algorithms are extremely sensitive towards the starting points as well as initial conditions [7]. The 
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implementation of stochastic behavior in evolving algorithms based on population helps in achieving an 
optimum solution. A survey is conducted on smart techniques application for voltage control in reactive 
power systems [8]. Scholars such as Kennedy and Eberhart worked on developing a progressive and 
powerful computation technique named PSO which is based on a simple concept and can be executed with 
few lines of coding. This computer code is based on primitive operators of mathematics and is cheaper to 
execute in terms of computational memory and speed [9]. This PSO technique is already in use for problem 
optimization in power systems and is also [10], [11], executed for resolving the popular issues on economic 
dispatch. For power optimization issues, Yoshida et al. [12] has worked on a two-stage method for executing 
a criterion for voltage stability. In the first stage of the method, the PSO helps in reducing the losses of 
transmission power. In the second stage of the method, the highest loading constraint of the system is gauged 
with help of continuous power flow method abbreviated as CPFLOW [13]. The initial solution fulfilling the 
basic requirements of system is identified as the ultimate best solution. Nevertheless, implementation of the 
PSO is limited for power market scheduling, particularly the reactive power market which is based on criteria 
for voltage stability [14]. 


2. PARTICLE SWARM OPTIMIZATION 

In the year 1995, Kennedy and Eberhart developed a search technique based on a population named 
as particle swarm optimization (PSO) method [15]. It is derived from the swarm’s natural movement in space 
searching for food. PSO can be used in optimizing mutidimensional problems, liner and nonliner problems 
and also single-dimensional problems. These algorithims are executed based on the particle population. 

The problem dimension is identified based on variables present in the problem. The particles in the 
search space are tagged with initial values. With an assumption of an initial velocity, these particles move in 
search space adjusting their position with respect to global best and personal best points in search space. In 
each iteration, the component of velocity changes, hence the initial value assumption is crucial [16]. Another 
important component is inertia weight. Finally, after tracking the complete search space, the particle reaches 
the best position. 

The parameters on which the PSO algorithm depends are inertia component, velocity factor and 
retardation factor. Hence it is crucial to identify the best value for these parameters [17]. To find the global 
optimum position, search space is explored extensively based on PSO parameters such as particle’s 
acceleration, velocity and inertia weight. If the values of these parameters are incorrectly assigned, the 
algorithm may not provide complete coverage and get trapped into a local optimum point. At time denoted as 
t, the i” particle X;(t) can be written as: 


Xi = [Xn (©). Xp (the ove eee eee eee ne eee Xp (t)] (1) 


2.1. Basic of PSO model 
Where X, and V, are optimized parameters, X;,(t) is position of the i" particle with respect to the k” 


iteration, V;x(t) is velocity of the i" particle with respect to the k” iteration as shown in Figure 1. The particle 
changes its velocity and positions according to the following [18]: 


yt! =W vf T Cp * rand () * (Piesti — pë) + Cy ` rand () = (Onesti E g”) (3) 


go = x Ea p (4) 





Figure 1. Search technique of the particle swarm optimization (PSO) 
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Where, C, is self-confidence range, C, is known as the swarm range, V;* is velocity of the i''particle at 
iteration ‘k’, k is iteration number, w is inertia constant, Pye; is value of personal best position of each 
particle in i" iteration, gbesti is value of personal best position of each particle in i" iteration and rand () is 
random number between 0 and 1. 


2.2. PSO parameters 
2.2.1. Size of swarm 

The size of Swarm is identified based on count of particles in population. This count is associated 
with the performance of PSO. If there are a very less number of particles, the algorithm is swifted and 
trapped at a point of local optimization whereas if there is a very large count of particles, algorithm is speed 
is reduced. There is no particular rule to achieve the idle number of particles but the swarm size is increased 
with an increase in problem dimension [19]. 


2.2.2. Number of iteration 

Based on the problem, the maximum number of iteration is identified in particle swarm 
optimization. If the iteration count is less, algorithm is stopped at a premature stage whereas if the iteration 
count is more the algorithm complexity is increased along with an increase in computational time which 
slows the convergence [20]. In this paper, number of iteration is considered as stopping criteria which is 500 
iterations. 


2.2.3. Components of velocity 

The position of the swarm is based on the component of velocity which allows the particle to 
randomly moved according to particle’s best position and global best points. To start with, a random velocity 
value is tagged to each particle which is then updated in each iteration. The velocity value can be updated 
based on (3) for each iteration. 


2.2.4. Coefficients of acceleration 

The trajectory path on which the particle moves is determined by the coefficients of acceleration. 
This is affected by the experience of particles and other particle in search space. For an efficient tracking of 
global optimal position and convergence of the complete solution, the coefficients of acceleration must be 
properly assigned [21]. 


2.2.5. Inertia weight 

Another important factor in PSO is inertia weight determines the search space exploration along 
with convergence. Interia weights are available in several types such as linearly decreasing, random 
decreasing, and logarithm. Intertia weight which is decreasing linearly helps in improving the performance of 
PSO. Several types of inertia weights are examined with the iteration count, characteristics of convergence, 
algorithm complexity for a comparative analysis [22]. 


3. COMPUTATIONAL PROCEDURES 

The block diagram for PSO is shown in Figure 2 on the PSO optimization of load flow for reactive 
power compensation [23], [24]. There are fourteen steps need to be performed to get the results for 6 bus Oman 
electrical system as shown in Figure 3. Table 1 and Table 2 show the line data and bus data respectively. 


Figure 2. PSO block diagram 
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Figure 3. 6-bus Oman electrical grid 


Table 1. Line data for a 6-bus Oman electrical grid 


Line No. Bus Number Length (km) Total Impedances Total Impedances Half line charging Tap 
Q/Length (km) p.u/ Length(km) Admittance p. u ratio 

From To R X R X B/2 - 

1 1 3 20 0.4283 2.82100 0.002458 0.01619 0 - 

2 3 4 15.5 0.3319 2.18620 0.001905 0.012547 0 - 

3 4 5 16.5 0.3533 2.32730 0.002027 0.013357 0 - 

4 5 6 60.9 1.30417 8.58994 0.007484 0.049299 0 - 

5 6 2 34.5 0.738817 4.86622 0.004240 0.027928 0 


Table 2. Line data for a 6-bus Oman electrical grid 


Bus No. Voltage Generation Load 
V (pu) O(deg) P(MW) Q(MVAR) PMW) Q(MVAR) 
1 1.00 0 - - - - 
2 1.00 0 70 0 - - 
3 1.0 - - - 115 38 
4 1.0 - - - 0 0 
=) 1.0 - - - 50 17 
6 1.0 - - - 102 34 


The steps PSO optimization of load flow for reactive power compensation are as clarified below: 
NBUS = number of buses and this paper are Oman 6 bus electrical grid. 

Calling Y bus function. 

Setting base MVA = 100. 

Calling bus data function as provided by Oman electricity transmission company. 
PSO parameter initilization. 

— Number of particles = 66. 

— Maximum number of iterations = 500. 

Creating storage matrix for variables. 

Initializing the variables. 

— Voltage magnitude = random values between 0.9 and 1.1. 

— Load angle = random values between 0.3 and -0.3. 

— Velocity of voltage = random values between 0.5 and -0.5. 

— Velocity of theta = random values between 0.5 and -0.5. 

h. Initializing the index for checking reactive power limit as in Table 3. 


ono of 


om rm 


Table 3. Index for checking reactive power limit 


SI.No. Condition PVIND 
1 Q calc < Q min 1 
2 Q calc > Q max 2 
3 Q min < Q calc < Q max 0 


Particle swarm optimisation for reactive power compensation on Oman 6 bus ... (Adnan Saif Al Mamari) 


1916 O ISSN: 2088-8694 


pà 0 


Computing the initial values for all particles. 
For i= 1 to maximum number of iteration, initialize inertia weight until reach the stopping criteria 
which is 500 iterations. 
Update the values of voltage and angle. 
Calling load flow calculation PROGRAME to calculate line loss and power injected at each bus. 
. Calculate slack bus power which bus 2 as shown in Figure 3. 
End. 


Qi © 


Pe BTM 


4. RESULTS AND DISCUSSION 

As explain above, the algorithm of PSO is implemented on 6-bus Oman electrical system in Matlab 
by using the bus and line data as in Table 1 and Table 2, and the result as discussed in Table 4 and Table 5 
below. The Newton Raphson (NR) as conventional method has been implemented on the same 6 bus 
electrical system as in [25], and the comparison between results is highlighted as well. 


Table 4. Load flow result of 6-bus Oman electrical system by using PSO 


Bus No. Voltage (P.u) Angle (Degree) Injection Generation Load 

MW Mvar MW Mvar MW Mvar 

1 1.000 0.000 479.550 -423.23 479.55 -423.23 0.00 0.00 

2 1.000 -2.287 -309.460 69.58 201.93 69.58 0.00 0.00 

3 1.060 -4.762 289.350 293.98 -194.46 331.98 115.00 38.00 

4 1.080 -5.950 -533.360 785.7 289.35 785.7 0.00 0.00 

5 1.000 -8.490 -533.360 -447.35 -483.36 -430.85 50.00 16.50 

6 0.973 -5.439 -104.470 -123.67 -19.47 -98.67 85.00 25.00 
Total 23.530 155.012 273.5 234.512 250.00 79.50 


Table 5. Line power flow result of 6-bus Oman electrical system by using PSO 


Bus Bus Line Loss 
From To P(MW) Q(Mvar) From To P(MW) Q(Mvar) MW Mvar 


1 3 479.55 -423.23 3 1 -469.49 489.46 10.06 66.23 

3 4 160.04 -195.48 4 3 -158.96 202.6 1.08 7.12 

4 5 448.3 583.1 5 4 -438.92 -521.25 9.39 61.85 

5 6 -94.44 73.9 6 5 95.52 -66.83 1.07 7.07 

6 2 -199.99 -56.84 2 6 201.93 69.58 1.93 12.74 
Total Loss 23.53 155.012 


It is obvious that, the real power loss by using PSO is 23.530 MW which is more than what is found 
by using the NR method which is 1.35 MW. The computational time for PSO (0.69 sec) is less than in NR 
(1.45 sec). In this implementation, the number of iterations is for PSO, in search for optimized reactive power 
compensation is 500 iterations and it is also used as the stopping criteria. In case of increasing this number, 
the computational time of convergence will increase as well. The parameters dependency in PSO algorithm is 
obtained using different inertia weight by changing different parameters such as tolerance factor, retardation 
factor, cognitive factors. The comparison between PSO and NR is shown in Table 6 and in terms of voltage 
magnitude is as in Figure 4. 


Table 6. Comparison between NR and PSO 
Criteria PSO NR 
Number of iterations 500 5 
Total power loss (MW) 23.53 1.35 
Computational time(sec) 0.698793 1.45 
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Figure 4. Comparison of voltage magnitude for NR method & 
PSO method 
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5. CONCLUSION 

In this paper, the PSO algorthim was used to obtain the power loss for 6-bus Oman electrical system 
which will help to determine the reactive power compensation needed to maintain the required voltage 
profile. In addition, the results have been compare with NR results for the same system and it is clear that 
PSO is more robust and effective. PSO is more flexible and has better and fast convergence compared to NR 
method and other techniques. 
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